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Abstract 


Surface albedo measurements have, in the past, been 
taken mainly near the earth's surface and from aircraft. 

In this study a surface albedo was determined from visible 
data collected by the NOAA-4 polar orbiting meteorological 
satellite. In order to filter out the major cause of 
atmospheric reflectivity, namely clouds, techniques were 
developed and applied to the data resulting in a map of 
global surface albedo. Neglecting spurious surface albedos 
for regions with persistent cloud cover, sun glint effects, 
insufficient reflected light and, at this time, some 
unresolved influences, the surface albedos retrieved from 
satellite data closely matched those of a global surface 
albedo map produced from surface and aircraft measurements 
and from characteristic albedos for land type and land use. 
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Introduction 


The albedo of an object is a measure of its light- 
reflecting property. Usually expressed in percent, the 
surface albedo is the ratio of the amount of reflected 
radiation per unit area to the total amount of incoming 
incident radiation per unit area. 

Reflectivity is a function of the nature of the 
surface, the angle of incidence, and the angle from which 
the measurement is taken. For liquids such as the oceans, 
reflectivity varies with depth, turbidity, surface roughness 
and current velocity. Albedos for solids such as snow, 
vegetation, and soil, vary with color, texture, wetness, 
grain size, and, in the case of snow, age. Therefore, due 
to the inhomogeneity of the surface, the earth's albedos 
range from small values, e.g., 6 to 9% for oceanic regions 
near the Equator, to large percentages, e.g., 70 to S5% for 
fresh fallen snow. Both changes in the solar zenith angle 
and (because of non-diffuse reflective surface properties) 
changes in the angle from which the reflective measurement 
is made also affect the measured albedo. 

Prior to the advent of meteorological satellites, 
albedo measurements were taken exclusively near the surface 

m 

and from aircraft. These produced both regional and global 
maps of surface albedo (Kung, Bryson and Lenschow, 1964; 
Posey and Clapp, 1964) . 
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With the introduction of satellites as meteorological 
observatories^ global maps were first produced from 
vidicon (television) camera measurements (Winston and 
Taylor# 1967; Taylor and Winston# 1968; Winston# 1971) and 
then from radiometric measurements (Raschke and Bandeen# 

1970; Vonder Haar and Suomi, 1971; Gruber# 1973; Flanders 
and Smith# 1975) . However# these maps are of planetary rather 
than surface albedo# the difference being that planetary 
albedo is a measure of the earth-atmospheric reflectivity 
which includes effects of reflectance, scattering and 
absorption by air molecules# aerosols and clouds. 

The major atmospheric influence on satellite measure- 
ments in the visible spectrum is the highly reflective 
nature of clouds. It is the aim of this study to develop 
techniques for filtering out cloud reflectance from 
satellite measurements in order to retrieve surface albedo 
quantities. Moreover# a global surface albedo array 
derived by these techniques has been constructed from five 
days of satellite data and evaluated in comparison to other 
such arrays and characteristic surface albedo values. These 
techniques have been computerized by the author so that 
weekly# monthly# seasonally, or yearly climatologies of 
global surface albedo can be generated. 

Development of surface albedo arrays would contribute 
to many areas in meteorological research. For instance# 
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the surface albedo Is one of the essential components 
in radiation and heat budget investigations. Since 
differential heating of the earth's surface provides the 
energy to drive the circulation of the atmosphere, the 
global surface albedo distribution is directly related to 
atmospheric circulation energetics. In atmospheric 
simulation models, such as the Goddard Institute for Space 
Studies (GISS) global atmospheric circulation model 
(Sommerville et al., 1974), a global surface albedo is cn- 
of the input arrays necessary to run the model. Better 
surface albedo specification would presumably lead cu 
better results in these areas. 

On the question of whether or not global atmospheric 

models are sensitive to changes in surface albedo, tests 

have been conducted at GISS with the GISS global 

atmospheric model by Drs. W. Quirk and Y. Sud (1976) and 

analyzed by the author. Two five day forecasts were 

produced using the same initial conditions, but different 

global surface albedo arrays, A climatological January 

global surface albedo array from Schutz and Gates (1972) 

was used in one forecast and a daily updated array in the 

other. This latter array was also initialized from the 

Schutz and Gates array but changed to correspond to the 

weekly average snow and ice boundary map produced by the 

1 

Satellite Applications Group of NESS (NOAA) . It was 

l.The National Environmental Satellite Service (NESS) of 
the National Oceanic and Atmospheric Administration. (NOAA) 


updated daily using a functional relationship between the 
model's forecast snow depth and surface albedo. After 
five days, surface temperature differences of 5* to 15 # F 
and sea level pressure differences of 4 mb, occurred in 
both the U.S, and Canada in regions of systematic, non-random 
differences in albedo. These results indicate not only 
model sensitivity to albedo variations, but also a possible 
impact on forecast performance. 

Data Source 

The data for this study are derived from measurements 
taken by the scanning radiometer aboard the NOAA 4 
meteorological satellite ( Environmental Satellite Imagery , 
Feb, 1976 ; Fortuna and Hambrick, 1974; Conlan, 1973; 

Schwalb, 1972). The NOAA 4 is a sun-synchronous, near 
circular, polar-orbiting satellite completing one orbit 
every 115 minutes. It transits the Equator on the 
descending node (southbound) at approximately 0900 hours 
local time and at approximately 2100 hours local time on 
the ascending node (northbound) . 

The two-channel scanning radiometer measures energy 
in both the visible (0.5 to 0.7 jam.) and the infrared 
(10,5 to 12,5 pm,) spectral ranges, scanning from horizon 
to horizon across the orbital path or track. Preflight 
calibration consisted of correlating the response of the 




radiometer to a brightness source of known spectral 
energy output* The only onboard calibration is a baseline 
determination when the radiometer scans empty spice. 

For the visible channel the instantaneous field of 
view is a 4 km. square at the nadir (subsatellite point) , 
increasing to a rectangle 7.5 km. by 15 km. at a horizontal 
distance 1668 km, from the nadir. At the subsatellite 
point there is a 4 km, gap between successive cross track 
swaths of data which disappears more than 1385 km, from 
the subpoint. Swaths of data from successive orbits overlap 
1668 km. from nadir (satellite observation angle of 
approximately 60°) at the Equator. This distance from the 
subpoint decreases with increasing latitude, Torward 
movement of the spacecraft combined with crosstrack 
scanning provides global visible measurements every 12 or 
13 orbits. 

The brightness measurements taken by the visible 
channel of the scanning radiometer are recorded by an 
onboard magnetic tape recorder. This analog signal is later 
transmitted to one of two ground receiving stations and 
retransmitted to NESS in Suitland, Maryland for processing. 

Processing includes earth location, normalization, 
digitation, cropping and mapping. Each data spot must be 
earth located in order to determine the solar zenith angle 
at the time of observation for that spot. The data are then 


normalized to an overhead sun using a cosine function 
factor of the solar zenith angle. The analog signal i3 
converted into digital brightness counts from zero to 254. 
This is a linear conversion whereby each digital count 
-epresents a range of 40 foot lamberts (i.e., digital 
brightness count 0 represents 0 to 39 foot lamberts, etc.). 
Digital brightness count 255 is used for missing data. 

Gaps in the data stream occur periodically off the 
Eastern U.S. coast near Bermuda due to insufficient data 
recording capabilities of the onboard recorder, the 
location of the two ground receiving stations, and shifting 
orbital tracks. 

For mapping purposes, one day of data, beginning 
with the orbital pass where the satellite transits the 
Equator at approximately 070 0Z, is grouped together. In 
regions where data are collected from more than one pass, 
the latest data are retained. This results in lines of 
discontinuty on visual displays between data retained 
from successive passes, and especially between data 
retained from the first and last pass of the day. The 
digitized data are then placed into two 2048 X 2048 
square grids, each overlaying a hemispheric polar 
stereographic map. These are made available on magnetic 
tape through the Satellite Data Services Branch of the 
National Climatic Center, Data for five days, January 31 


through February 4, 1976, have been used in this study. 

At the time this data set was being processed, the 
digital code for missing data, 255, was not being inserted 
into polar grid points which had a solar zenith angle equal 
to or greater than 90°, Any brightness recorded by the 
radiometer, whether real or spurious was left intact in 
the data stream. This resulted in non zero digital counts, 
usually multiples of 10, for grid points in the polar night 
region. Photographic and facsimile maps exhibit u spotty 
appearance in these regions due to this non zero field. 

This same spotty appearance occurs on maps over 
portions of Antarctica and the Weddell Sea, northern 
Russia and Scandinavia for the days studied. Again the 
reason for this is an insufficient amount of reflected 
light. These areas are scanned on the last orbital pass 
for that day, thus all the data collected, no matter how 
far away from the satellite (i.e., large satellite 
observation angles) are mapped. Any previous measurements 
for this region from earlier passes are simply deleted. 

Since the satellite overflies a region in the morning its 
crosstrack scanner measures regions "toward" the sun and 
"away" from the sun. Regions viewed "away" from the sun 
would have a greater solar zenith angle. This larger solar 
zenith angle coupled with the retention of all data on the 
last pass leads to saving data collected in the polar regions 
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with insufficient illumination. 

Method of Analysis 

The data on tape are in digitized counts of reflected 
brightness from both the surface and the atmosphere. To 
retrieve a surface reflected brightness the atmospheric 
reflected brightness must be filtered out. Since clouds 
are the m«.jor contributors to atmospheric reflection of 
visible light , techniques were developed and applied to the 
data to filter out this cloud effect. No attempt was made 
at this time to filter out the less influential 
atmospheric effects such as reflectivity due to atmospheric 
aerosols or Rayleigh scattering. 

The earth's surface was assumed to be an isotropic 
surface (i.e., with no preferential direction of reflected 
energy) thereby eliminating the effect of measurements 
taken at different satellite observation angles and solar 
zenith angles. This assumption was also used by Winston and 
Taylor (1967), Taylor and Winston (1968), and Flanders and 
Smith (1975), for studies of brightness, long wave 
radiation and albedo. 

The cloud filtering technique was based on the premise 
that clouds increase the amount of light reflected back 
to space. For example, over a dark surface the brightness 
count for a particular grid point is higher on a cloudy 
day than on a cloudless day. In this study, groups of grid 



point brightness counts were analyzed in the form of 
histograms* An increase in the modal brightness value 
is generally observed on cloudy versus cloudless clays. 

As sky conditions change from clear to thin cirrontratus 
to lower, thicker stratus, the modal brightness also shifts 
to higher brightness values; Overcast stratiform cloud 
conditions over the analyzed region gave only one . 
brightness peak* Open and closed cellular cumulus, 
cumulonimbus and regions only partially overlain by clouds 
resulted in bimodal histograms. In this later ^ase. the 
lower end peak, (i.e., the peak with the lower brightness 
count) and adjacent brightness counts were taken to 
represent the surface reflectivity. The higher end peak 
presumably represented cloud reflection. 

The first cloud filtering technique employed was 

developed to eliminate any higher brightness peaks and save 

only the lower surface representative peak. In order to 

carry out a low end modal searching process the hemispheric 

data (2048 X 2048 grid points) were regrouped into quare 

grid boxes consisting of 80 grid points on a side 

(approximately 800 km.). This formed a 25 X 25 grid box 

array for each hemisphere. A histogram of brightness counts 

2 

versus frequency was then constructed from the 6400 grid 

point values for each grid box. As an example, figure 1 

is a histogram constructed from one grid box over the 

2, A special computer program was designed to automate 
construction of the histogram, as well as the search 
for modal brightness, and determination of characteristic 
surface albedos. The Fortran code for the program is give 
in the Appendix, 
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. Fig.l. Histogram of brightness counts versus smoothed frequencies 
for one grid box over the Indian Ocean. The lowest, middle.' and 
highest row of numbers are the brightness counts, actual frequencies 
and smoothed frequencies respectively. Each star represents two 
smoothed frequencies. 


Indian Ocean. The bottom row of numbers are the brightness 
counts from zero to 255. Above these are the corresponding 
frequencies or number of values of that brightness measured 
within the grid box. For instance# for brightness value 
16 there are 264 measurements of a 16 brightness count out 
of 6400 total measurements. A three point smoothing operator 
of the following form was used on brightness values 2 
through 253 to reduce the "noise" in the histogram: 

SF n = 1/3 (F n »i + F n + F n+ D 

where 

SF n = smoothed frequency for brightness count n. 

F n = frequency of brightness count n. 

F n+ 1 = frequency of brightness count n + 1. 

Again using brightness count 16 as an example# the 
frequencies of brightness count 15, 16 and 17 (369# 264 and 
258 respectively) are added together and then divided by 3 
to give the smoothed frequency for brightness count 16 
(i.e.# 297). A two point smoother was applied to brightness 
values 1 and 254 with the next higher and lower brightness 
counts respectively. 

In other words the actual frequencies of brightness 
counts 1 and 254 were added to brightness counts 2 and 253 
and then divided to 2 to give smoothed frequencies for 
brightness counts 1 and 254, 

The smoothed frequencies in figure 1 appear in a row 


just above the actual frequencies. The stars are used to 
give a more graphic representation of the smoothed 
frequencies for each brightness count. Each ctar represents 
2 frequencies up to a possible 240 frequencies (i.e., 120 
stars) . 

Initially a mode was determined from the lowest 35% of 
the data (i.e., frequencies being added together starting 
with brightness count 1 and reaching 35% of the data) • 
However, for regions of clear skies the low end mode did not 
correspond to the mode for the entire histogram. It did 
correspond to the lower end mode in the bimodal case. Modes 
from clear sky regions were studied and an empirical 
relationship was derived between the mode and the beginning 
of the peak (basically, the first brightness count with a 
frequency of 4 or more) . This range was then used as a 
criterion to distinguish modes of clear skies from modes 
of partly cloudy skies for modal values from more than 
35% of the data. Modes from the first 80% of the data and 
then from 10% less data down to 50% and then 5% less data 
to 40% were checked against the range criterion for clear 
skies. The first mode which had a range within the limiting 
criterion or the mode for the lowest 35% of the data was 
saved. This mode was then used in calculating a mean 
brightness. 

In figure 1 a low end modal search process was carried 


out for the first 80% of the smoothed frequencies. The 
modal brightness value was 21 with a mode of 330 smoothed 
frequencies. The range between the modal brightness value 
and the first brightness value of four or more (i.e., 
brightness count 1 in this case) is 21 which is less than 
the clear skies criterion and so no further low end modal 
searching with successively smaller percentages of data was 
conducted. Next the number of brightness counts between 
the mode and the beginning of the brightness peak was 
added to the modal brightness count and used to define the 
end of the brightness peak. In figure 1 the number of 
brightness counts between the modal brightness count of 
21 and the beginning of the brightness peak, namely 
brightness count 1, were computed. This range of 20 was 
then a<!fed onto the modal brightness count in order to 
designate the end of the brightness peak. Brightness count 
41 was thus defined as the end of the brightness peak in 
this case. A mean brightness count was calculated from 
smoothed frequencies of brightness count 1 to the brightness 
count marking the end of the brightness peak. The mean 
brightness count was 21.8 for figure 1. These mean 
brightness counts represent the mean surface brightness for 
that grid box. Surface representative means were calculated 
for each grid box for each of the five days. 

These means, however, are in terms of brightness counts 
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representing specified internals of energy expressed in 

foot lamberts. To convert these nteans into albedo 

measurements , each value must be multiplied by a constant 

conversion factor derived by Gruber (1974), It takes into 

-1 

account a solar constant of 1.94 ly.min. and the percentage 
of total incoming solar radiation effectively sensed in 
the visible spectral range of the radiometer. 

The first cloud filtering technique was used to capture 
a surface representative mode for clear or partly cloudy 
sky conditions* A second filtering technique had to be 
applied to the data to eliminate means from regions with 
overcast conditions where the satellite could not effectively 
"see" the earth's surface. A five day minimum scanning 
procedure based on the work of McClain and Baker (1969) in 
the reduction of cloud contamination in snow brightness 
measurements and snow boundary delineation for satellite 
data was used. A similar multiday minimum scanning 
procedure was used by Chen (1975) for monthly minimum 
planetary albedo mapping from satellite data. 

From five days of calculated mean albedos for each grid 
box only the lowest mean albedo was retained as the surface 
representative albedo. In this final scanning procedure 
the means had to have been calculated from a sample of at 
least 2000 frequencies , that is, the number of smoothed 
frequencies from brightness count 1 to the brightness count 
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marking the end of the modal peak had to exceed 1999. This 
eliminated the cases where the low end mode occurred very 
close to the beginning of the peak thus giving a lower than 
representative mean. 

Results 

Hemispheric surface albedo maps generated from the 
application of two cloud filtering techniques on five 
consecutive days, January 31 to February 4, 1967, are 
given in figures 2 and 4. The 25 X 25 grid overlaying the 
geographical hemispheric map shows the boundaries of each 
grid box. The lowest surface representative mean albedo 
(which passes the limiting criteria) for the five days is 
written in each box. An I for any box denotes insufficient 
light occurred in the northern hemisphere in the polar night 
region and in both hemispheres due to the retention of all 
data from the last pass. A 99 found in four of the grid 
boxes in the northern hemisphere means that no surface 
representative albedo was determined from the data. The 
reasons for this will be discussed later. 

For most grid boxes the albedos are close to 
characteristic surface albedos for specific types of surfaces. 
Albedos for water surfaces generally fell in the 5 to 12% 
range with higher albedos in the higher northern latitudes. 

Non snow and ice covered land regions had albedos of from 
7 up to 35% in desert regions. Snow and ice covered areas 
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exhibited higher albedos than its characteristic non 
covered form* Albedos for Antarctica were from 80 to 
91% * It must be noted that the scanning radiometer was 
only able to distinguish a reflective brightness up to a 
corresponding 91% albedo. At this point the instrument was 
"saturated". 

Besides a relative closeness of most computed albedos 
to the characteristic surface albedos , some computed values 
showed discrepancies. For these grid boxes daily hemispheric 
photographs produced from the same scanning radiometer data 
were studied. The first problem was one of persistant cloud 
cover over the grid box for the five days. This did not 
necessarily mean an overcast condition over the entire grid 
box* Boxes with as low as 7 tenths sky cover produced 
unusually high albedos. These albedos have been labeled 
on the maps with a bar to denote persistant cloud cover and 
do not represent a surface albedo. 

Other grid boxes with unusually high albedos were a 
result of clouds and sun glint. Sun glint occurs only yver 
the water and appears as a bright area. This mirror j,ifect 
happens when the viewing angle is equal to the angle of 
incidence of the sunlight. In the photographs this takes 
the form of bright swaths along the satellite track or 
bright spots. Albedos recorded for these regions are 
generally around 35%. A combination of cloudy days, missing 
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data for one or more days, and the occurrence of sun glint 
on the day with the lowest mean albedo for that region, 
produced higher than expected albedos. 

The reason for unusually higher albedos from some grid 
boxes is not as clear. In these cases clouds were present 
covering less than 7 tenths of the grid box area on the day 
of lowest mean albedo. These are denoted by a filled in 

I 

triangle in the lower right corner of the grid square. 
Further study needs to be conducted to determine why the 
high albedos are present and how to correct this situation 
in the computer program. 

One further point should be kept in mind when viewing 

the surface representative albedo maps. When a grid box 

covers two or more types of surfaces the lower albedo 

(provided enough cases for that surface are present) 

will represent the entire grid box. The most numerous 

example of this is for grid boxes covering both land and 

water. If the water has a lower albedo as is usually the 

0 

case then the albedo for the water would represent the 
entire grid box. However, if clouds covered the water for 
the five days then the land albedo would be retained, 

A comparison between global surface albedo arrays 
derived from characteristic albedos for surface types, 
from readings taken near the ground and from aircraft, 
and the global surface albedo array derived from satellite 
data follows. 
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Posey and Clapp (1964) compiled four monthly average 
global surface albedo mans drawing exclusively upon the 
work of Budyko (1958) for open ocean values and defining 
characteristic albedos for different land use types and 
surface conditions for non open ocean areas. The average 
January surface albedo map, figures 3 and 5 will be used 
for comparison. 

Southern Hemisphere 

Discounting grid boxes with persistent cloudiness, 
sun glint and unresolved errors, the satellite albedos 
derived for the open oceans exhibit a fairly good agreement 
with Posey and Clapp. The range is from 5 to 12% as opposed 
to 6 to 9% for Posey and Clapp, The satellite derived maps 
show more geographical fluctuation owing to the fact that 
these are five day minimums rather than monthly means. The 
systematic stepwise increase in albedo found in Posey and 
Clapp is not present in the satellite derived maps. 

Antarctic albedos are generally higher, 80»to 91% 
than the near uniform 80% albedos of Posey and Clapp. Four 
grid boxes in the eastern section of the continent which 
cover both land and water depict the case where the land 
albedo is higher than the water since the water is covered 
by persistent clouds. The retrieved albedo would then be 
the lower, cloud albedo for these cases. Thus these boxes 
are marked for further study. 



Fig. 2. Surface albedos derived from satellite data, 
January 31 - February 4, 1976, Southern Hemisphere. 
Albedos are given in percent. A bar over an albedo 
indicates persistent cloud cover over five days. A 
double asterisk above an albedo indicates sun glint 
contamination. A * in the lower right corner of the 
grid box means that the albedo was unrepresentative 
for reasons as yet' unclear. 


































Fig. 3. Normal surface albedo, January, Southern Hemisphere. 
Solid lines in Antarctic waters are northern boundary of 
pack ice (heavy) or boundary of ice regions (light). No 
snow' boundary indicated.- Numbers arc aJ.bedo in percent. 
(Taken from Posey and Clapp) . 
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Australian albedo 3 do not show the large range of 
albedos found in Posey and Clapp. numerous reasons for 
these differences such as grid size, solar zenith angle at 
time of observation, uncharacteristic soil wetness, changes 
in land use since 19C4, etc., could be ascribed for these 
differences but the cause (s) is unclear at this tine. 

Agreement is very close, 7 to 13% compared to 7 to 10%, 
for most albedos in Southern Africa. The portion of South** 
western Africa with albedos from 18 to 30% is not apnearant 
on the satellite derived map. 

Slightly higher albedos, up to 5% occur, over most of 
South America except along the western coast. 

Northern Hemisphere 

Most open water albedos are again close to those given 
by Posey and Clapp. Fluctuations of albedos for the same 
latitude still exist, however, at high latitudes the increase 
of albedo with increasing latitude is found in both maps in 
the North Atlantic and North Pacific oceans. 

Uigh albedos from 49 to 59% due to sea ico are found in 
three grid boxes: one covering the Bering Straits, another 
extending from the southwestern edge of Greenland across to 
Labrador and a third in the Shelikhov Gulf, north of the 
Sea of Okhotsk between Siberia and the Kamchatka 
Penninsula » In order to verify the existence of sea ice 
in these areas at the time the observations were made, 



Fig. 4. Surface albedos derived front satellite data, 
January 31 - February 4, 1976, Northern Hemisphere. 
Albedos are given in percent. A bar over an albedo 
indicates persistent cloud cover over five days. A 
double asterisk above an albedo indicates sun glint 
contamination. A m in the lower right corner of the 
grid box means that the albedo was unrepresentative 
for reasons as yet unclear. An albedo of 99 represents 
missing data. 
















‘Fig. 5. Normal sur.’ .oe albedo, January, Northern Hemisphere. 
Numbers arc albedos in percent. Heavy dashed line over 
northern oceans is boundary separating ice navigable 
to ice breakers from that navigable to heavy ships: over 
land, snow boundary (or line of 50 probability of snow) . 
Light solid lino mainly in Artie- Ocean, boundary of solid 
pack ice. Light solid lines over land are discontinuities 
in albedo, following land-use boundaries. 

(Taken from Posey and Clapn). 


weekly average snow and ice boundary (WASIB) maps for the 
periods January 26 to February 1, 1976 and February 2 - 8, 
1976 are included, figures 6 and 7. These maps are produced 
by the Satellite Applications Group of NESS (NOAA) using 
photographs derived from data from NOAA-4 and SMS-1 
satellites. Sea ice boundary extents are denoted by circles. 
The lowest representative surface albedo for the Bering 
Straits grid box was recorded on January 31, 1976 when sea 
ice was present. Ice conditions changed the following week 
according to figure 7. Ice conditions for the other regions 
remained the same for the two week period. 

Four 99 values appear in the Northern Hemisphere. A 
value of 99 means that no surface representative albedo 
was recorded for any of the five days. This resulted 
when an insufficient number of cases comprising the low end 
peak was present for each of the five days. Mostly cloudy 
skies over the five davo woulu cause such a condition. 

Surface representative albedc • for Central America, the 
northern portions of South America and Africa are in close 
agreement with those of Posey and Clapp, The higher albedos 
for the Sahara desert stand out well in both maps, with the 
satellite derived values slightly higher, 35% as opposed 
to 30%. 

The albedos are also in close agreement in the Arabian 
Penninsula except for grid points covering both land and 
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February 2-0, 1976. 




water* For these the lower albedos for water predominate. 

Albedos are somewhat lower in Iran, Afghanistan and 
Pakistan for the satellite derived maps and slightly 
higher in India, Southeast Asia and the non desert areas of 
China. The desert regions of China, exhibit albedos of 
20 to 30% in Posey and Clapp whereas they range from 16 to 
24% in the satellite derived albedo maps. 

The snow boundary extending from the Caspian Sea to 
the Pacific Ocean on the WASIB maps and on Posey and Clapp, 
denoted by a dashed line, lie close to each other except 
the former displays snow cover in the Himalayan Mountains. 
Albedos to the north of the snow boundary are similar with 
40 and 50% albedos near the boundary and decreasing to the 
upper 20' s, 30' s and 40% albedos farther to the north on 
both maps. 

In Europe the snow extent on the WASIB maps is more 
widespread than on the Posey and Clapp map. Higher 
albedos of 31 to 47% occurred in the non snow covered 
regions of 13% in Posey and Clapp. 

The snow boundary line for North America was fairly 
close for the Posey and Clapp and the WASIB maps. Judging 
from the depth of mow on the ground maps, for February 2 
and 9, 1976 (figures 8 and 9) extensive melting of snow in 
the midwest and western plateau regions was indicated for the 
week preceding February 2, 1976, with widespread snowfall in 



Fig. 8. Depth of snow on the ground, February 2, 1976. 

(Taken from Weekly Weather and Crop Bulletin, February. 3, 1976). 



Fig. 9. Depth of snow on the ground, February 9, 1976. 

(Taken from Weekly Weather and Crop Bulletin, February. 10, 1976). 
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the following week. Melting snow and snow free areas 
characteristically have lower albedos than the same regions 
with fresh fallen snow. As expected albedos in the north 
central U.5, and Rocky Mountain states have lower albedos , 

16 and 17%, than 45 and 55% for Posey and Clapp. The non 
snow covered southern states have identical albedos. The 
grid boxes in the southwestern states and south into Mexico 
display differences possibly linked to ground wetness and 
a mixture of soil types. Except for Hudson Bay, both maps 
exhibit similar albedos over Canada, Alaska's albedo of 
38% is the same for Posey and Clapp except near the southern 
coast where Posey and Clapp note a 62% albedo region. 

Conclusions 

The cloud filtering techniques applied to satellite 
scanning radiometer data in the visible spectrum enabled 
the production of a global surface albedo map with quite 
reasonable values. Open water albedos were generally 7 to 
11% and increased to 21% in the higher northern latitudes. 

Three regions with extensive sea ice had albedos of 49 to 
59%, Non snow and ice covered land albedos ranged from 
7 to 16% with albedos from the mid 20s to 35% for desert 
regions. Snow covered lands displayed high albedos from 
the mid 20s to 63%. Albedos for Antarctica were normally 
high ranging from 80 to the radiometers upper limit of 91%. 

Incorrect albedos did, however, still remain in the 
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global albedo map. Disregarding regions of insufficient 
reflected light due to the tilt of- the earth and data 
retention procedures used on the last pass of the satellite, 
some unusually high albedos were caused by persistent clouds 
over the regions for five days, by a combination of clouds 
and sun glint, and by other not as clearly defined influences. 
Further study needs to be undertaken using more days for 
minimum surface albedo retrieval, and a smaller minimum 
sample size than the 2000 used in this study for a 
representative surface albedo determination for each grid box. 

Nevertheless, a comparison of global albedo maps, one 
derived from satellite data and the other from surface and 
aircraft measurements and characteristic albedos for specific 
land use and land types showed two quite similar maps for 
all surface conditions. 
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GO TC E3 

165 XMfc ANsSCM/NUM 

AM CANS ( XVE*N *0. 2 60 I ) 

0IFF-NCy*CLV-SUM**2 
C 1 V = n I F F ✓ ( MM* ( MJM— 1 ) ) 

FDLV=5LRI l Cl V) 

ASDE VI - ( XJ- E AN- SCE V)*C.36CI 
A SCti V-(XVL AN+ SCE V ) *C. 36 C I 

V.MI 11(6.1 /C) AVfcAN. ASOE VI . ASDE V2 • I H I SB • MOCEBV . IHI SE . MCOESF.NUM, 

l rci_M 

1 7 G FDRMAK/. AX** NEA N AL EE C t = * .. Ft» • 2. 4X . • AL EEOC RANGE ( 1 ST AN.UEV) =• . . 

1 rt.i. ' TC • .F 7.2/. AX . • IH I SB= ♦ .14 .* MCO E ev= • . I 4 . • IHISE=*.I4, 

2 * V.on&SF= • . 1 5. * M.<V=*.I5.« P CENT CF OAT A SEARCHED** «F4 .2/J 

WR I Tfc ( <J ) I F STDA . AWE AN » < IAFRECtLH.NI ).LH = l ,25 6 ) . 

1 t 1 SFEEQtL 1 .N l > .LI = l .25 0 , ASDE VI . AS0EV2. 

2 NUM . : H l S3 . IHISE.MCOEU V , NCDESF .PCENT 
53 COM TIM. H 

V.R nt [ 6. PC JNLKREC .lOACC.RA.l'V/.Le.IVALE.I IV ALB. I I VALE 
BO FCRPATC/./.* LAST RECCPD RE A C= * . IS . 3X . • I 640 0=* . I £ , 3X . • NAa » , X 4, ✓ . 
t • I V.- LH- • . I 5. 3 X • * 1 V AL t = • ,15 *JX . *1 I VALB= *. 15. 3X.*IIVALt-= •*!£./. 

25 CONTINUE 
ENDK ILL 9 

PAUSE * MCLNT NEXT TAPE(FF) CN LOG 8* 

IF 5 TCA= IF S IDA* 1 

.FI I F ST D A • LE • 1 1 S3 ) GO TC 6 

STOP 

END 


